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Abstract 

We present an analysis of twist-2, leading order QCD amplitudes for hard exclusive 
leptoproduction of mesons in terms of double/nonforward parton distribution func- 
tions. After reviewing some general features of nonforward nucleon matrix elements 
of twist-2 QCD string operators, we propose a phenomenological model for quark 
and gluon nonforward distribution functions. The corresponding QCD evolution 
equations are solved in the leading logarithmic approximation for flavor nonsinglet 
distributions. We derive explicit expressions for hard exclusive ir°, 77, and neutral 
vector meson production amplitudes and discuss general features of the correspond- 
ing cross sections. 
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1 Introduction 



Our current knowledge of the sub-structure of nucleons is based to a large extent on high- 
energy scattering experiments which probe its quark and gluon distribution functions. 
The most prominent processes in this respect are deep-inelastic scattering and Drell- 
Yan leptoproduction. In addition a large amount of information can be deduced from 
measurements of electromagnetic and weak form factors. 

Recently new observables, namely nonforward parton distribution functions, attract 
a great amount of interest. Although being discussed already some time ago [I], [!|, 
they were introduced in the context of the spin structure of nucleons recently in Ref . |4| . 
Nonforward parton distributions are a generalization of ordinary parton distributions 
@, H . However they are also closely related to nucleon form factors. Thus they combine 
different aspects of the nucleon structure and offer new insights. 

To summarize the relation of nonforward parton distributions to ordinary parton dis- 
tributions recall that at the twist-2 level the latter can be represented as normalized 
Fourier transforms of forward nucleon matrix elements of non-local QCD operators ||, 
constructed as gauge-invariant overlap of two quark or gluon fields separated by a light- 
like distance. Nonforward parton distributions are defined by the same non-local QCD 
operators - just sandwiched between nucleon states with different momenta and even- 
tually spin. Equally close is the relation to nucleon form factors which are defined by 
nonforward matrix elements of the same QCD operators taken in the local limit. 

Nonforward parton distributions are probed in processes where the nucleon target 
recoils elastically. To select twist-2 correlations a large scale has to be involved. One 
possible process is deeply- virtual Compton scattering as investigated in Ref.||. Another 
promising class of reactions sensitive to nonforward distributions is hard exclusive lepto- 
production of mesons (for recent works see |5|, [7|, |J and references therein). This is based 
on a factorization theorem which has been proven recently in Ref. H. It concerns the kine- 
matic domain of large photon virtualities Q 2 ^> Aq CD and moderate momentum transfers 
\t\ ~ Aq CD . The theorem states that for incident longitudinally polarized photons the 
meson production amplitudes can be factorized in a perturbatively calculable part, de- 
scribing the interaction of the virtual photon with quarks and gluons of the target, and 
matrix elements which contain all information about the long-distance non-perturbative 
strong interaction dynamics in the produced meson and the nucleon target. The latter 
are nothing else than the wanted nonforward distribution functions. 

Data on the exclusive production of neutral vector mesons have been taken lately 
at high center of mass energies at HERA (for a review and references see |9[]). In the 
measured kinematic domain the corresponding production cross sections are controlled 
by the nonforward gluon distribution of the target. The latter has been approximated in 
earlier investigations by the ordinary gluon distribution of the target nucleon, [[K], [11], |12|, 
O. The quality of this procedure is still a matter of debate (see e.g. Ref.|TH||). 

At smaller center of mass energies, typical for current and future fixed target experi- 
ments at TJNAF ESI, HERMES EE and COMPASS [T7L a similarity of the nonforward 
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and ordinary gluon distribution in the accessible kinematic domain is not likely. Further- 
more also the nonforward quark distribution of the target will contribute in a significant 
way. Therefore information about nucleon nonforward matrix elements can be obtained 
from these experiments. In addition, measurements of other exclusive processes, like 
pseudoscalar or charged vector mesons production, will allow to variety of new 

nonforward quark and gluon distribution functions ||. 

In this paper we outline the calculation of hard exclusive leptoproduction amplitudes 
for a variety of mesons. To obtain first insights in the magnitude and behavior of the 
corresponding cross sections a modeling of nonforward parton distributions is necessary. 
We use a simple ansatz which satisfies constraints due to forward parton distributions, 
form factors and discrete symmetries. Furthermore, we discuss the Q 2 -dependence of non- 
forward quark distributions in the nonsinglet case by solving the corresponding evolution 
equation. We are then in the position to provide baseline estimates for the differential 
cross sections for 7r°, 77 and neutral vector meson production. 

The structure of this paper is as follows: in Sec. 2 we review the definition and gen- 
eral features of double and nonforward distribution functions. A simple model for double 
distributions is presented in Sec. 3. In Sec. 4 the Q 2 -evolution of nonsinglet quark distri- 
butions is outlined. Exclusive meson production amplitudes are derived in Sec. 5. After 
presenting results in Sec. 6, we summarize in Sec. 7. 



2 Double and nonforward distribution functions 

We first recall definitions and basic features of double and nonforward distribution func- 
tions. For this purpose we stay close to the notation of Ref. 0. In leading twist quark 
and gluon distribution functions of nucleons are defined by forward matrix elements of 
non-local light-cone operators sandwiched between nucleon states ||: 

(N(P, S)\ 6(0, z) \N(P, S)) z2=0 ~ fdx [e-^f(x) ± e^/(x)] . (1) 

J 

Here \N(P, S)) represents a nucleon with momentum P and spin S. The operator O(0, z) 
stands for a product of two quark or gluon fields, respectively, separated by a light-like 
distance z ~ n, with n- a = a + = a + a 3 for an arbitrary vector a. In Eq.([TJ) we suppress 
any dependence on the normalization scale. Note that the sign in front of the anti-parton 
distribution / is determined by the charge conjugation property of the operator O in 
the local limit. Forward matrix elements as in Eq.fll]) are explored in processes where 
no momentum is transfered to the target, e.g. in forward virtual Compton scattering or 
deep-inelastic scattering, respectively. 



2.1 Double distribution functions 

In exclusive production processes a non- vanishing momentum r = P — P' is transfered to 
the target. Thus in leading twist rather nonforward matrix elements of quark and gluon 
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light-cone operators are probed. These can be parametrized through so called "double 
distribution functions" which where introduced in Ref.||. For the unpolarized quark 
distribution one has: 

(N(P', S') | ^(0)S [0, z] ^(z) \N(P, S)) z2=0 = 

N(P', S') z N(P, S) f dx f X dy \ e -*t(P-*)-*u(r-*) > F ( Xj y ^ t ) _ gi^P-*)-^™)^ ^ ^ 



+N(P', S') a ^ff U N(P, S) Cdx [ X dy\e- ix(p - z) - iv{r - z) K(x,y,t) - e^^^ K(x, y, t) 
iM Jo Jo 1 

(2) 

where we use the notation z = z^ . The operator on the LHS of Eq.([|) is built from 
quark fields connected by the path-ordered exponential [0, z\ = Vexp[—igz,j, Jq 1 d\A' M (zX)} 
which guarantees gauge invariance and reduces to one in axial gauge n ■ A = (g stands 
for the strong coupling constant and denotes the gluon field). On the RHS we have the 
product of the Dirac spinors for the initial and final nucleon. The quark and antiquark 
double distribution functions F, F, K and K depend on the momentum transfer t = r 2 , 
and two light-cone variables x and y (x — 1 — x, y — 1 — y). Their interpretation is 
discussed at length in Ref . |J . Note that the distributions K and K enter proportional to 
the momentum transfer r. 

A comparison of Eqs.(|l|) and (0) demonstrates that in the forward limit, r — > 0, the 
double distributions reduce to ordinary quark distributions, e.g.: 



f(x) = / dyF(x,y,t = 0). (3) 
Jo 

For the polarized quark and the unpolarized gluon distribution one has |J: 



(N(P\ S')\m^ 5 [0, z] $(z) \N(P, S)) z2=Q - 



N(P', S')z>y 5 N{P, S) dx dy e- ix(p - z) - iy{r - z) AF{x, y, t) + e" (p ' 2) -^ (r ' 2) AF(x, y, t) 
Jo Jo L 

+K — terms, 

(4) 

z fl z v {N(P',S')\TrG^(0)[A;0}G^(z),\N(P,S)) z2=o = 

N(P', S')zN(P, S)^-^ I dx [ X dy\e- ix(p - z) - iy{r - z) + e^-*)-^*-*)] G(x,y,t) 

4 jo Jo 1 J 
+K — terms, (5) 



with P = (P' + P)/2. In Eqs.([|,[Sp we have only indicated the so-called "i^-terms" which 
are proportional to the momentum transfer r as in Eq.([|). Of course, relations similar to 
(|3p hold also for the double distributions in Eqs.(f|,[5]). 

It is important to realize that the double distributions defined above fulfill a symmetry 
constraint based on hermiticity. For the unpolarized quark distributions F and F this 
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can be seen by considering the C-even (E) and C-odd (A) combinations: 



E 
A 



} 3 ^ 



z z 
2'~2 



if* 



(6) 



From the definition (Q), combined with the fact that double distributions are real, it 
follows that E and A are purely imaginary and real, respectively: 



E = 2iN(P',S')zN(P,S) 
-i 



x dx dy F(x, y, t) + F(x, y, t) 
Jo Jo 1 

A = 2N{P',S')zN{P,S) 



x dx dy F(x, y, t) - F(x, y, t) 



sin 



cos 



xP ■ z + r ■ z \ y — 
xP ■ z + r ■ z [y — 



x 



2 

1 — x 



+ K — terms, 

+ K — terms. 

(7) 

As a consequence E = (E — E^) /2 as well as A = (A + At) /2 is symmetric with respect to 
an exchange of the initial and final nucleon momentum, i.e. P <-» P' and simultaneously 
r <-> — r. This yields the following sum rules: 



dx j dy F(x, y, t) + F(x, y, t) cos (xP ■ z^j sin 



dx / dy F(x, y, t) — F{x, y, t) 



sin ( xP ■ z ) sin 



r • z yy 
r ■ z [y ■ 



1 — x 
2 

1 — x 



0, 
0. 



(9) 



For the polarized quark nonforward distributions as well as for the gluon distribution 
similar constraints can be derived. They imply that the double distributions in Eqs. (H,§,0) 
are symmetric with respect to an exchange of 



1 — x — y. 



(10) 



This symmetry, besides being important for modeling double distribution functions, is 
crucial for establishing proper analytical properties of meson production amplitudes. 

In the forward limit, r — >• 0, Eqs. (|2|,|4]||) immediately reduce to the definitions of 
ordinary twist-two parton distributions ([I]). On the other hand in the limit z — > 
Eqs.([|,f|,[5D define familiar nucleon form factors. One therefore obtains the following sum 
rules for double parton distributions R f|, EN: 



dx 



dx 



dy [F(x,y,t)-F{x,y,t)\ = F 1 (t), 
dy \k(x, y, t) - K(x, y, t)] = F 2 (t), 



dx dy AF(x,y,t)+AF(x,y,t) =G A (t), 
Jo 1 J 

dx / dy G(x,y,t) + G(x,y,t) = Q(t). 



(11) 
(12) 

(13) 

(14) 
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Fx, F 2 and Ga are the Dirac, Pauli and axial- vector form factors of the nucleon. The 
gluon form factor Q(t) is experimentally not observable. Its dependence on the momen- 
tum transfer has been estimated within QCD sum rules fllEfl , and is determined by a 
characteristic radius of the order of 0.4 — 0.5 fm. Similar sum rules hold for the remaining 
distributions [01. 



2.2 Nonforward distribution functions 

In the definitions of double distribution functions (Hdll) the variable y always enters in 
the combination x + y ( = X with ( = j^. As a consequence it is possible to define e.g., 
unpolarized nonforward quark and antiquark distribution functions ||: 

F c (X,t)= f 1 dx f dyF{x,y,t)5{X-{x + C,y)) 
Jo Jo 

= 6(X>C) [ X/C dyF{X-y(,y,t) + Q(X<0 f dy F{X - y(,y,t), 
Jo Jo 



F c (X,t)= dx dyF(x,y,t)5(X-(x + (y)) 
Jo Jo 

= 6(X>C) [ X/C dyF(X-y(,y,t) + e(X<() dy F(X - y(,y,t), (15) 
Jo Jo 

and similar for the polarized quark and gluon case. Since the scattered nucleon is on its 
mass shell, i.e. P' 2 = M 2 with the nucleon mass M, one finds £ < 1. Together with the 
kinematic constraint x+y < 1 this gives < X < 1. The variable X can be identified with 
the nucleon light-cone momentum fraction of the parton being removed by the operator O 
from the target, while the light-cone momentum fraction of the returning parton is equal 
to X - c. 

In terms of the nonforward distribution functions the light-cone correlators in (^|^J^) 
read: 



(N(P',S')\iP(0)H0,z}^(z)\N(P,S)) 



z 2 =0~ 



N(P', S') z N(P, S) f dX \e- lX{p - z) F c (X, t) - e l(x ~ c)(p - 2) F C (X, t) 



+N(P', S 



e -ix(P-*) K( ( X ,t) - e^ x -^ p ^K c (X,t) 



(N(P>,S')\mzi 5 [0,zMz)\N{P,S)) 



z 2 =0~ 



N(P', S')zj 5 N(P, S) f dX \e~ iX ^ AF C (X, t) + e i(x -° (p - z) AF ( (X, t)] + K - terms 



(N(P', S')\ TrG^(0)[A; 0]G e »(z), \N(P, S)) z2=0 - 



N(P\ S')zN(p,s)^[ 1 dX 
4 jo 



e -iX(P-z) + e i(X-0(P-z) 



G c (X,t) +K- terms. (16) 
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Note that in the literature also a second parametrization of nonforward distributions is 
popular [|J. It is based on a different asymmetry parameter £: 

( = iTrz = 2=? with os?s1 ' (17) 

Using for example for the unpolarized quark and antiquark distributions 

F'(u,£,t) = ! dx [* F{x,y,t)6{u-{x + £{x + 2y-l))), 
Jo Jo 

F'{u,£,t) = I dx f X F(x,y,t)5(u+(x + £(x + 2y-l))), (18) 
Jo Jo 

allows to define an alternative nonforward distribution as: 

H(u,$,t) = e(u>OF\u,^t) + e(Z>u>-0(F'(u,Z,t)-F'(u,Z,t)) 

- e(-£>u)F'(u,Z,t). (19) 

The parametrization of the nonforward matrix element (|2|) in terms of the new distribu- 
tions reads [§J: 

(n(p', s')\mz to, z) m \n(p, s)) z2=0 = 

N(P', S') z N{P, S)e~ i£ ? J du H{u, f , t)e~ iu{p - z) +K- terms. (20) 

Both definitions of nonforward parton distributions are of course equivalent. In the fol- 
lowing we will always parameterize nonforward matrix elements as in Eqs. ([15| , [T6|) , i.e. in 
terms of variables X and (. 

3 Model 

Nonforward parton distributions have not yet been measured. Therefore one has to rely 



on models []19j in order to provide estimates for exclusive production cross sections. To 
guarantee a proper analytic behavior of the involved amplitudes, as given by dispersion 
relations [20|, it is favorable to model double distributions (@,§||) instead of nonforward 



distributions. We parametrize the former, denoted generically by F(x, y, t; /Xq), such that 
they fulfill all constraints we are aware of. As a model ansatz at some low normalization 



scale /1q we take [21 



F(x,y,t;i4) = ij^hix,y)f(t), (21) 

where f(x, [Aq) stands for the corresponding ordinary quark and gluon distribution, respec- 
tively. It is clear that one can model in this way only the distributions F, F, AF, AF and 
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G. The "if-terms" in remain unconstrained because the corresponding forward 

parton distributions are not known. We choose: 

h(x,y) = 6y(l-x-y), (22) 

such that F (x , y , t\ Hq) satisfies the symmetry constraint in Eq. flHX) . Of course there are 

r -i 2 

many possible choices for h(x,y), e.g. h(x,y) = 12 y — . The latter however results 
in a nonforward parton distribution which is not continuously differentiable at X — £. 
Although we are not aware of any principle which forbids such a behavior, we use in the 
following the parametrization from Eq.(|22|) leading to nonforward distributions smooth 
in X. 

The form factor f(t) in Eq. (|21|) is responsible for the t-dependence of the double 
distribution functions. Motivated by the relationship between double distributions and 
nucleon form factors in Eqs.(|ll| - |TJ|) we assume: 

'»=Mtf) 2 ' (23) 

For quark distributions we take A from fits to the nucleon vector and axial form factors 
2"2|j . It should however be mentioned that the meson production amplitudes, as derived 



in Sec. 5, are controlled by nonforward distribution functions with different C-parity as 
compared to the above nucleon form factors. In the case of the gluon distribution the 
scale A can be taken from the QCD sum rule estimate in Ref. ||18||. One should keep in 
mind that a factorization of the t-dependence into a form factor f(t) is reasonable only 
at small values of \t\. At large \t\ the minimal component of the nucleon light-cone wave 
function dominates p3fl . This should lead to a change of the entire x and y dependence of 
F which cannot be absorbed in a pre-factor. However, as long as we are interested in the 
small- \t\ region, the parametrization should be acceptable from a phenomenological 
point of view. 

In Fig.|Tj we present our phenomenological model for the nonforward polarized w-quark 
distribution, Au^(X, //q), taken at t — 0. Here and in the following we use parametriza- 
tions for the corresponding forward distributions, which enter through Eq. (pl|) , from 
Ref. ||24j| . One can recognize the characteristic distribution amplitude like shape for X < £, 



and the similarity to ordinary parton distributions at X > 



4 QCD evolution of nonforward distributions 

An important ingredient for the discussion of nonforward distribution functions is their 
evolution, or dependence on the normalization scale. In QCD it is governed by renor- 
malization group equations which are direct generalizations of the well known DGLAP 
equations for ordinary parton distributions. 

Instead of looking at the distribution functions themselves it is convenient to consider 



QCD evolution at the operator level 0, E5[. Evolution equations for various quanti 



ties, like evolution equations for meson distribution amplitudes [23, 26], or the DGLAP 
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equations for ordinary parton distributions, follow then by taking appropriate operator 
matrix elements. At the one loop level evolution equations for operators can be solved by 
an expansion in a set of multiplicatively renormalizable operators which are determined 
through the conformal symmetry of one-loop QCD [pjj. 

In the following we outline a convenient numerical method to perform the QCD evolu- 
tion of nonforward parton distributions. Details of the method and a thorough discussion 
of its accuracy will be given in a separate publication. In the nonsinglet case the Gegen- 
bauer moments of nonforward parton distributions 



C f (n;tfj) = [ 1 dXC 3 J 2 (2X/C~l)^(X;^ ) 

•J 

are multiplicatively renormalizable ||, and their evolution reads: 

In 



(24) 



C c (n;/i ) 



C c (n;/x ), 



(25) 



with (3 = 11 - |iV/, The nonsinglet anomalous dimensions 7 n are given as usual by: 



(n + l)(n + 2) 



ra+l i 
j=2 J 



(26) 



Because the Gegenbauer polynomials C^j 2 form an orthogonal set on the segment [—1,1], 
the straightforward inversion of Eq.(|24]) is possible only for £ = 1. On the other hand one 
can expand the Mellin moments of the nonforward parton distribution 



f 1 dXX n FAX;nl) 
Jo 



(27) 



in terms of Gegenbauer moments Q(n;/Xo), and evolve the latter according to Eq.(|25| 
The corresponding expression has been obtained in Ref. f|: 



M n (C;/i 2 



n (—l) k 



M fc (C;/io) 



x 



E 



(-l)'2(2/ + 3)(fc + / + 2)! 
(n + l + 3)\(n-l)\(l-k)\ 



2d 



:28i 



To obtain the evolved distribution from its moments one usually performs an inverse 
Mellin transformation. In the present case, however, Eq. (|28|) has a simple form only for 
integer values of n, and therefore usual contour integration methods (see e.g. ||28|| ) in 
the complex n-plane are difficult to implement. This problem can be circumvented by 
expanding the nonforward parton distribution F^(X; /i 2 ) in terms of e.g., shifted Legendre 
polynomials Pk(2X — 1) [29] which are orthogonal on the segment [0, 1]: 



F c (X;/i 2 



£(2A; + 1) c k (Cfi 2 )P k (2X 

k=0 



(29) 
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The coefficients Ck{(; fJ 2 ) can be determined from the set of moments M n ((; /x 2 ): if a 
shifted Legendre polynomial Pk{2X — 1) follows an expansion 

P fe (2X-l)=5>£X", (30) 

n 

the coefficients Ck((; /i 2 ) are given by: 

c fc (C;^ 2 ) = E< M -(C;^ 2 )- (31) 

n 

In Fig.^| we present typical results for the evolution of the nonforward polarized 
li-valence quark distribution, parametrized at /j% = 4 GeV 2 according to our model 
in Sec. 3. For the leading order evolution we use Nf — 4 and Aqcd = 250 MeV. 
The evolved distribution behaves as expected: the area under the curve remains con- 
stant as a consequence of the vanishing anomalous dimension 70. The evolution shifts 
the distribution to smaller values of X, approaching very slowly the asymptotic shape 
An c (X;/i 2 ^oo) = 6C(0)X/C 2 (l-X/C) §. 1 



5 Exclusive meson production amplitudes 

A solid QCD description of hard exclusive meson production is based on the factorization 
of long- and short-distance dynamics which has been proven in Ref . M . The main result is 
that the amplitudes for the production of mesons from longitudinally polarized photons 
can be split into three parts: the perturbatively calculable hard photon-parton interac- 
tion arises from short distances, while the long-distance dynamics can be factorized in 
terms of nonperturbative meson distribution amplitudes and nucleon nonforward parton 
distributions. In the following we derive production amplitudes for a variety of neutral 
mesons in terms of these building blocks, staying at leading order in the strong coupling 
constant a s . 

The kinematics is defined as follows: the initial and scattered nucleon carries a mo- 
mentum P and P' = P — r, respectively, while t = r 2 stands for the squared momentum 
transfer. The momentum of the exchanged virtual photon is denoted as usual by q with 
Q 2 = —q 2 . Then q' = q + r is the momentum of the produced meson. Calculating the 
hard subprocess to leading twist accuracy allows to neglect the momentum transfer t, 
the invariant mass of the nucleon target P 2 = P' 2 , and the mass of the produced vector 
meson q 2 as compared with the virtuality of the photon Q 2 . Then the Bjorken scal- 
ing variable XBj = Q 2 /2P ■ q coincides with the longitudinal momentum transfer, i.e. 

( = r ■ n/ P ■ n = XBj- 

In leading order in the electromagnetic interaction we obtain for the S-matrix in ques- 
tion: 

S = (M^Nt mS^Uq) N(P)), 

1 Recently a similar method has been proposed independently in Ref. [j30[ . 
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= -i Jd*xe- iq - x (M(q')N(P')\e L - J(x)\N(P)). (32) 



Here |7£(?) N(P)) denotes the incoming photon and nucleon while {M(q') N(P')\ repre- 
sents the final meson and nucleon. The polarization vector of the longitudinally polarized 
virtual photon is given by e£ = ^(q' + XBjPY, and is the electromagnetic current. 
In second order perturbation theory we arrive at: 

S = 27cia s Jd 4 x d*y d 4 z e~ iq - x (M(q')N(P') \M(x, y, z)\N(P)), (33) 

with the time-ordered product of quark and gluon fields: 

M(x, y,z)=T (f(,) 7 ^ fc ( Z )f(x)4f(x)^( 2 /) 7 V(|/)i^) A?(yj) t A ah t B de , (34) 

where t A and t B denote the generators of color SU(3). As a next step we perform the usual 
Operator Product Expansion of A4(x,y, z). Choosing, for example, the flow of the hard 
momentum q as in Fig.^|(a), we obtain a term which can be interpreted as photon- meson 
transition in a background quark field provided by the nucleon. In the remaining (not 
shown) leading order contributions the photon couples to one of the other possible quark 
lines. A different choice of the hard momentum flow results in a photon-meson transition 
in the background gluon field as shown in Fig.|^(b). This gluon contribution has been 
discussed several times in recent papers [7|, [31]. We therefore outline in the following 
the calculation of the quark part. Here the diagram in Fig^(a) corresponds in leading 
order a, to: 



M(x,y,z) 3 \i> d (y)Yr(y)] \r{x)e L S^ a \x-z)Y^\z)] Sff A \y 



(35) 



S q and S 9 stand for the quark and gluon perturbative Feynman propagators, respectively. 
Next we perform a Fierz transformation in color- and Dirac-space. Of course in color- 
space only the singlet piece survives when matrix elements between colorless hadronic 
initial and final states are taken. In Dirac-space we have: 



# d (v) r r(y)} [r(x)e L s q(ca) (x - Z ) y ^\z) 

with the operators Oi and Oj defined as: 



52C°fd i ( y ,z)d j (x,y), 



(36) 



Oi,j{y,z) = ip(y)Tijip(z) 



(37) 



The coefficients C« are given by: 



16 



x 



^)7 p r, 



with r„-e{l,7 5 ,7 M ,757 M ,^}- 



(38) 
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For color singlet final states higher order corrections generate path-ordered exponentials 
which ensure gauge-invariance of the operators Oij ||. As a next step leading twist pieces 
have to be extracted. In the massless limit, i.e. P 2 = P' 2 = q 2 = 0, this can be achieved 
by using the Sudakov decomposition of 7 matrices: 

Y = 1 f 1 P" + 1 ^q ni + ^- (39) 

P-q P ■ q 

Typically only longitudinal terms proportional to P or q' contribute at the twist- 2 level. 
Evaluating the corresponding trace and sandwiching operators Oj and Oj between initial 
and final states gives the desired amplitudes in terms of meson distribution amplitudes 
and nucleon nonforward parton distributions. An important outcome of this calculation 
is that only operators with Tj = Tj contribute. Since the spin and parity of the produced 
meson determines the Dirac matrix in the meson amplitude, it automatically fixes the 
involved nonforward distribution. 



5.1 Vector mesons 

To leading twist accuracy two amplitudes for longitudinally and transversely polarized 



vector mesons exist (see [32"| for a recent discussion). In the longitudinal case one has: 

(M(q')\$(x)PiP(y) |0> = q' ■ P fy L C dr^y L {r) C *-(W), (40) 

Jo 

where fy L and §y L denote the corresponding decay constant and distribution amplitude. 
According to our previous discussion this implies that the quark contribution to the 
production of longitudinally polarized vector mesons is determined by the nucleon matrix 
element (N(P')\tjj(z)q'ifj(y)\N(P)) . As a consequence the nonforward quark distributions 
F and K, defined in Eqs.(PJl6|), enter. The corresponding amplitude reads: 



A q 



Ira* 



C F 1 N{P', S')q'N(P, S) 



N C Q 



P-q' 



Jo 



(It 



TT 



X 



dX 



F c (X,t)+F c (X,t) 



1 



+ 



X-ie X 
C F 1 N(P',S')(q'r-rq')N(P,S) 



x 



dX 



IMP -q 1 

K c (X,t)+K c (X,tj 



1 



-h 



+ 



C + ieJ 

1 



TT 



X-ie X-C + ie 



(41) 



Since the virtual photon and the final state vector meson carry the same C-parity, 
vector meson production can take place also in the gluon background field. Its contribution 
is |, 0,0: 



A 9 

Ay L 



2n a* 



1 1 N{P', S')q'N(P, S) 



N C Q 



P-q' 



( ) f 

TT JVl 



TT 



11 



1 dX (Y ^F^VXM + K ~ t6rmS ' ^ 
o [X — ie){X - Q + le) 

with G being the nonforward gluon distribution (||) of the nucleon. 

Transversely polarized vector mesons are described by the amplitude 



(M(q')\^(x)a^(y) |0> = -iM - ^Jfe f dr <S> Vt (t) {rx+fy) . (43) 



o 

Their production from longitudinally polarized photons involves chiral-odd nucleon non- 
forward parton distributions. Thus, in general, only the nonforward quark transversity 
distribution contributes. An explicit calculation shows, however, that in leading order 
a s the twist-2 contribution vanishes and one has to go either to higher orders, or to 
higher-twists to obtain a non-zero result. 

5.2 Pseudoscalar mesons 

In the case of pseudoscalar mesons one ends up with the amplitude: 

(M(g')l ^(x) 75 A%) |0> = iq'-Pf PS f 1 dr $ P5 (r) e i( ^ + ^, (44) 



with the meson decay constant fps and the distribution amplitude $>ps- According to 
the above discussion the production process is sensitive to the nonforward nucleon matrix 
element (N(P')\ifj(z) , -f5q'if}(y)\N(P)) . Thus the nonforward polarized quark distributions 
AF and AK from Eqs.(f|,|T6D enter. Collecting all leading order contributions gives for 
the pseudoscalar meson production amplitude: 

, • C F l N(P',S'h 5 q'N(P,S) j; yi , $ P5 (r) 
^P5 = z7ra s — — — fp S dr 



x 



N C Q P-q> J ^ Jo rf 

f 1 dX \AF C (X, t) - AF ( (X, t)} (^— + v * ) +K- terms, (45) 
Jo L J \ A — it X — C + te I 



Due to C-parity conservation pseudoscalar meson production receives contributions only 
from C-odd nonforward distribution functions. As an immediate consequence two gluon 
exchange is impossible. At least three gluons have to be exchanged, which in the language 
of twist expansion corresponds to higher- twist. 



6 Results 

Using the derived production amplitudes together with the model distributions from Sec. 3 
allows to calculate the production of neutral pseudoscalar and vector mesons. In both 



cases we use the asymptotic meson distribution amplitude [|3|, |26|1 : 

$P 5 (r) = $ yL (r) = 6r(l-r). (46) 

In numerical calculations we have neglected the practically unconstrained "iT-terms". As 
their contribution enters proportional to r it is bound to be small at small momentum 
transfers. 
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6.1 7r and 77 production 

To obtain predictions for a specific process the generic amplitudes from Sec. 5 have to be 
furnished with flavor charges. Using standard SU(3) wave functions for the pseudoscalar 
meson octet implies for ir° and rj production^ an replacement of the nonforward distri- 
butions in Eq . fl45|) through distributions with specific flavor as listed in Tab.l. For 7r° 





AF ( (X,t) - AFs(X,t) 






§ (Au ( - Aw c ) + | (Ad c - Ad ( ) 




V 




J (Au c - Au c ) - i (Ad ( - Ad ( ) + | (As c - As f )] 



Table 1: Flavor structure for 7r° and rj production. Analogous relations hold for the 
"K - terms". 



production we use the standard value for the decay constant f w = 133 MeV. In Fig.|] we 
present the corresponding differential production cross section for a proton target taken 
at t — t min = — x 2 B jM 2 /{I — %Bj)- We restrict ourselves to the region \t min \ < 1 GeV 2 , 
which implies XBj < 0.6. 

One finds that, up to logarithmic corrections, the production cross section drops as 
l/Q 6 - Furthermore, at small values of XBj the calculated cross section is proportional to 
Xgj with A ~ 0.6. This behavior is controlled by the small-xsj behavior of the polarized 
valence quark distributions which enter in Eq. fl2T|) . For the used parametrizations from 
Ref.p^l one has indeed XBjAu v (xBj) ~ XBjAd v (xBj) ~ x B j. The decrease of the pro- 



duction cross section at large xsj is due to the form factor in Eq.(21). It is important 
to note that this behavior should be quite general and largely independent of a specific 
model for double distributions. The reason is that the rise of the production cross section 
with increasing Xsj is brought to an end through the decrease of the involved nonforward 
or double distributions at large momentum transfers \t\. Since the latter is determined 
by a typical nucleon scale, A ~ 1 GeV, the maximum of the pseudoscalar meson pro- 
duction cross section at t = t min should occur when t min starts to become sizeable, say 
—t min /A 2 > 0.2. This implies XBj ~ 0.3 in accordance with our result. 

Our prediction for the 77 production cross section turns out to be approximately a 
factor 2/3 smaller than for 7r°, but of similar shape. This is due to the comparable shape 



of the polarized u and d valence quark distributions from Ref.|24| which enter in Eq. (|2T 



6.2 Vector meson production 



For p°, $ and u production the nonforward distributions in Eq. (f4l| , f42"D have to be modified 
according to Tab. 2 In the kinematic domain of small Xsj < 0.01 the quark part of the 

2 We restrict our considerations to the pure octet state rjs, neglecting mixing with the singlet state tjq. 
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F c (X,t)+F c (X,t) 


G c (X,t) 


p° 




j K + u c) + 1 K + d c)_ 






UJ 








172 G < 


$ 







Table 2: Flavor structure for p°, lo and $ meson production. 



production amplitude ([II]) is negligible. Then the well known SU(3) relation for the cross 
section ratios a(p) : <j(uj) : cr($) =9:1:2 follows immediately from Tab. 2. 

In Ref.[|i~0|, |i~T| , p~2[] hard exclusive neutral vector meson production has been investi- 
gated at small values of XBj within the framework of perturbative QCD. In this work the 
nonforward gluon distribution which, in general, enters the production amplitude ( |42| ) 
has been approximated by the ordinary gluon distribution. In the limit of small Xsj the 
obtained amplitude corresponds to Eq. ([42]) after the replacement G^((,t) — > (9(0: where 
g denotes the ordinary gluon distribution. However, note that for a typical behavior of 
the gluon double distribution G(x, y, t) ~ x~ x with A > one obtains from Eq.([3|) for the 

ratl °' Qc(C,* = o) i -h a Jo d yyG(Cy^ = Q ) > l m 

(9(0 ~ Jo dyG((,y,t = 0) ' 1 1 

For our model distribution from Sec. 3 this leads to an increase of the vector meson pro- 



duction cross sections by nearly a factor of two as compared to the results of [12 



It should be mentioned that the leading order calculations from Refs.[10, 11, 12] over- 
estimate the experimental data from HERA P] by a typical factor 5 — 10. Additional 
Fermi motion corrections have been added to fit experimental data. This emphasizes the 
need for a systematic investigation of higher twist effects |3~T| ]. 

In Fig.[5]we present the result for p° production from a proton in the kinematic domain 
of HERA as obtained from the amplitudes in Eqs.([|]],[|^), combined with the double 
distributions from Sec. 3. For the p meson decay constant f p = 216 MeV has been used. 
After multiplying with a suppression factor T(Q 2 = 16.9 GeV 2 ) « 0.1 from qualitative 
agreement can be achieved. 

Finally we would like to note that special attention should be paid to exclusive charged 
meson production. Here at large Q 2 nonforward parton distributions are probed which 
are related to matrix elements of twist-2 non-local operators off diagonal in flavor: 

d qq ,(x,y) = q(x)Tq'(y), (48) 

with flavor dependent quark fields, q = u,d, s. As a consequence hard exclusive charged 
meson production offers possibilities to explore exotic parton distributions which have 
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not been measured elsewhere. These obey isospin symmetry relations analogous to the 
Bjorken sum rule known from polarized deep- inelastic scattering: 

(n\O du (x,y)\p) = (p\O uu (x,y)\p) - (p\d M (x,y)\p), 

(n\6 du (x,y)\p) = (n\6 dd (x,y)\n) - (n\6 uu (x,y)\n). (49) 

Here (n\ and (p\ denote neutron and proton states, respectively. These sum rules imply 
a close relation of charged meson production cross sections. They allow, for example, to 
relate the isovector part of the amplitude for p° production to a linear combination of p + 
and p~ production amplitudes. 

7 Summary 

Ordinary parton distributions accessible e.g., in deep-inelastic scattering measure the 
nucleon response to a process where one parton is removed and subsequently inserted back 
into the target along a light-like distance, without changing its longitudinal momentum. 
Generalized parton distributions, so-called nonforward distributions, can be studied in 
deeply virtual Compton scattering and hard exclusive leptoproduction of mesons. They 
describe a situation where the removed parton changes its longitudinal momentum before 
returning to the nucleon. Furthermore, hard exclusive charged meson production provides 
possibilities to investigate new processes where the removed quark carries different flavor 
than the returning one. 

In this paper we have studied general properties of double and nonforward parton 
distributions. We have found a symmetry of quark and gluon double distribution functions 
based on their relation to nonforward matrix elements of QCD string operators. Its 
implication for phenomenological models has been outlined. 

Although the leading order QCD evolution equations seem to be more complicated in 
the non-forward than in the forward case, one can solve them explicitly using an expansion 
in a set of orthogonal polynomials instead of performing a Mellin integral in the complex 
plane of distribution function moments. We have outlined this method for the flavor 
nonsinglet case. A generalization to flavor singlet distributions is possible. 

In the second half of this paper we have derived the amplitudes for the hard exclusive 
production of neutral mesons from longitudinally polarized photons. After suggesting 
a phenomenological model for double or nonforward distribution functions which obey 
appropriate symmetries, we have presented results for exclusive tt°, r\ and p° production. 

In the pseudoscalar case nonforward polarized valence quark distributions enter. Sev- 
eral features of the presented results should be independent of our specific model for 
nonforward distribution functions: at small Xsj pseudoscalar meson production cross sec- 
tions drop for decreasing xbj- This is related to a similar behavior of the corresponding 
forward distributions. Furthermore, the production cross sections peak at moderate XBj- 
This is due to the dependence of the involved nonforward distributions on the momentum 
transfer which should be controlled by a typical nucleon scale. 
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We also have presented results for p° production in the kinematic domain of recent 
HERA measurements. Here the production process is dominated by contributions from 
the nonforward gluon distribution. For a large class of model distributions this tends to 
be larger than the corresponding forward distribution. 
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Figure 1: Polarized w-valence distribution An^(X;/i 2 = 4) from Eq.(|2~l~D for different 
values of (. 
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Figure 2: Evolution of the nonforward polarized w-valence distribution Au^o^iX; At 2 ) 
different values of /i 2 . 
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Figure 5: Cross section for p° production from a proton through the interaction of lon- 
gitudinally polarized photons. The data are taken from ZEUS |33| (open squares) and 
NMC H (filled circles). 
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